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INTRODUCTION 


This  is  the  final  report  that  describes  the  results  of  a  research  effort  into  the  development  of  nanosensors 
for  breast  cancer  biomarkers. 

THE  PROPOSAL  HAD  THREE  SPECIFIC  AIMS: 

1.  Design  of  a  cantilever  based  sensor  for  biomolecular  recognition,  using  Finite-Element 
simulation. 

2.  Fabrication  of  Biofunctionalized  Nanoscale  Sensors  capable  of  detecting  targeted  molecules  at  a 
concentration  of  less  than  1  ng/ml 

3.  Characterization  of  Functionalized  nanosensors  for  selected  breast  cancer  markers,  and 
comparing  with  existing  immunohistochemical  and  Fluorescence  in  situ  hybridization  (FISH)  techniques  on 
well  established  biomarkers  such  as  Her-l/neu,  estrogen  and  progesterone  hormone  receptors  in  tumor  tissue, 
and  selected  mucin  antigens  in  blood. 


We  report  that  we  have  achieved  in  full  two  of  the  specific  aims,  and  have  accomplished  part  of  Aim  3. 
We  have  characterized  the  funcationalized  nanosensors  and  have  detected  the  breast  cancer  biomarker  CA15.3 
in  solution  at  clinically  relevant  concentrations.  We  were  not  able  to  perform  studies  on  serum  samples  from 
breast  cancer  patients  for  the  reasons  mentioned  in  the  report  below.  Among  the  accomplishments  of  this 
project,  we  have  discovered  an  exciting  gating  principle  using  electrical  methods  that  has  the  potential  for 
significantly  enhancing  the  sensitivity  to  biomarkers  and  control  of  the  nanoscale  structures.  We  have 
published  this  in  several  papers  (References  1-4),  and  filed  a  patent  application  (Ref  5).  A  startup  company, 
Ninth-Sense  Inc,  is  negotiating  to  license  this  idea  from  our  University  with  the  goal  of  raising  fund  to  expand 
studies,  specifically  for  breast  cancer.  We  are  pleased  to  report  that  this  DOD  funded  project  has  resulted  in 
both  research  papers,  and  intellectual  property  with  the  potential  to  benefit  breast  cancer  research  and 
development. 


Key  Research  Accomplishments 

We  list  here  the  accomplishments: 

Task  1:  We  have  completed  simulations  to  characterize  and  develop  nanoscale  cantilever  sensors. 
Simulations  conducted  in  the  presence  of  water  show  that  hydrodynamic  effects  have  to  be  taken  into  account  in 
order  to  properly  characterize  the  sensitivity  of  these  sensors.  At  the  same  time,  we  have  performed  electrical 
conductance  measurements  on  nanonwire  channels  in  order  to  assess  the  sensitivity.  A  new  description  for  the 
electrical  response  of  the  nanochannels,  based  on  a  modified  Field-Effect  Transistor  model,  that  explains  the 
gating  principle  and  the  regime  of  operation  of  the  sensors  that  is  optimized  for  biosensing  for  breast  cancer 
biomarkers. 

Task  2:  We  have  successfully  used  electron  beam  lithography  to  fabricate  nanosensors  and  have 
functionalized  them  using  silanization  protocols.  Some  of  the  nanomechanical  structures  our  collaboration  has 
developed  are  among  the  smallest  and  highest  frequency  nanomechanical  resonator  ever  constructed.  In 
performing  conductance  measurements  on  the  structures,  we  have  discovered  a  very  interesting  gating  effect 
that  can  be  used  to  enhance  the  sensitivity  of  the  nanosensors. 

Task  3:  We  have  successfully  shown  that  the  nanosensors  respond  with  sensitivity  and  specificity  to  the 
breast  cancer  biomarker  CA15.3,  at  clinically  relevant  concentrations.  A  new,  more  sensitive  operating  regime 
for  the  semiconductor  sensors  has  been  discovered.  This  work  has  led  to  numerous  publications,  and  a  patent 
application  has  been  filed.  A  start-up  company  is  negotiating  for  the  license  to  commercialize  this  research  for 
breast  cancer  studies. 
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Reportable  Outcomes 
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Figure  1 :  (a)  Schematic  of  Silicon  nanowire  with  side  gates  and  electrodes.  The  nanowire 
is  exposed  on  three  sides  along  the  longitudinal  directions,  (b)  The  nanowire  shown  here 
is  300  nm  x  230  nm  x  8pm.  (c)  A  silicon  nanowire  with  an  Au/Ti  side  gate,  (d)  The 
scanning  electron  micrograph  displays  three  silicon  nanowire  devices  on  the  same  chip,  (e) 
An  optical  micrograph  shows  the  flow  chamber  sealed  on  top  of  the  devices  on  the 
interface  board,  (f)  Schematic  of  the  differential  conductance  (dl/dV)  measurement  circuit. 


We  have  used  simulations  to  characterize  and  develop  several  types  of  nanosensors.  Figure  1  (a)  shows 
a  schematic  of  the  nanochannel  sensor  developed  in  this  research  project.  Figure  1  (b-d)  show  electron 
microscope  images  of  a  single 
silicon  nanowire  with  contacts 
on  either  end,  and  more 
complex  devices  with 
fabricated  side  gates  for 
amplification  and  control.  Note 
that  the  width  of  the  sensors  is 
100  nm.  The  nanowire 
conductance  is  modulated  as 
selected  biomarker  targets  bind 
to  the  surface  of  the  nano  wire. 

The  steps  involved  include 
electron-beam  lithography, 
development  of  the  exposed 
pattern,  metallization  followed 
by  liftoff,  reactive  ion  etching, 
and  finally  metal  removal  to 
fabricate  the  device.  These 
procedures  are  standard 
methods  of  fabrication 

developed  in  the  semiconductor  industry  and  condensed  matter  physics, 
and  have  the  potential  to  be  mass  produced.  All  the  steps  have  been 
carried  out  in  our  laboratory  in  the  Physics  department  and  the  Photonics 
Center  at  Boston  University.  The  sensitivity  of  the  devices  arises  from 
the  large  surface-to-volume  ratio  at  the  nanoscale. 

The  surface  functionalization  protocols  are  based  on  silane 
chemistry,  and  have  the  ability  to  bind  to  any  selected  antibody  for 
specificity.  We  have  shown  that  an  important  new  step  of  Atomic  Layer 
Deposition  (ALD)  improves  the  efficiency  of  silanization.  Our  research 
work  has  shown  that  both  sensitivity  and  specificity  to  breast  cancer 
biomarkers  has  been  achieved  by  silicon  nanowire  sensors. 

Figure  2  shows  that  we  have  successfully  fabricated  an  array  of 
nanosensors  with  controlled  dimensions  for  highly  parallel  studies  on 
functionalized  nanowires.  These  images  demonstrate  that  lithography 
methods  have  the  potential  for  producing  special  purpose  biomarker 
chips.  Unlike  fluorescence  or  optical  methods,  the  chips  make 
measurements  purely  electrically.  The  semiconductor  industry  has  long 
shown  that  it  can  make  electrical  devices  smaller,  faster  and  cheaper.  In 
Figure  3(a),  an  optical  micrograph  picture  alsoshows  the  mating  to  a 
plastic  fluidic  channel.  The  channel  has  a  volume  of  about  35  pi,  and 
allows  for  fluid  exchange,  mixing  and  sequential  testing  using  a  syringe 
pump.  Unlike  other  microfluidic  channel  devices,  where  mixing  is  a  problem,  this  device  allows  mixing  to 
occur  in  the  sample  volume  just  above  the  nanowire.  The  electrodes  connecting  to  the  external  electronics  are 
insulated  from  the  serum  or  buffer.  The  entire  device  is  mounted  onto  a  semiconductor  chip  carrier,  for  easy 
integration  with  electronics. 

Figure  3  shows  the  device  in  operation  and  the  sensitivity  to  the  surface  charge  density.  Electrical 
characterization  of  the  device  showed  a  novel  operating  regime  of  reverse  bias.  This  new  operating  regime  was 
used  to  demonstrate  sensitivity  and  signal  amplification.  The  electrical  characterization  and  electron  beam 


Figure  2.  (Top)  Optical 
micrograph  of  the  device,  showing 
the  flow  chamber.  (Bottom) 
Electron  microscope  image  of  an 
array  of  nanosensors  of  different 
widths. 
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writing  was  performed  in  the  Physics  department,  and  the  rest  of  the  fabrication  was  carried  out  at  the  Photonics 
Center. 


TVtb  (trtn)  voltage 

Figure  3.  Electrical  characterization  and  demonstration  of  ultrasensitive  pH  detection  and  signal  amplification. 
(Left)  Differential  conductance  of  the  silicon  nanowire.  (Middle)  Change  in  the  differential  conductance  near  zero 
bias  as  a  function  of  pH  change  in  the  solution  in  the  flow  chamber.  Zero-bias  differential  conductance  is  used  to 
avoid  electrolysis  of  the  solution.  (Right)  Amplification  of  the  signal  by  the  application  of  local  gate  voltage  on  the 
nanofabricated  gate  close  to  the  nanowire  (see  Fig.  la  and  1c). 
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Device  configuration 


picture  of  device  from  top  view 


AI.O.  (22nm  thick) 
grown  by  atomic  layer 
deposition  (ALD) 


Al  as  top  gate 
Si  wire  ~100nm*100nm 
SiO2(380nm) 


Si  substrate  (~600pm) 


Using  these  devices,  we  have  been  able  to  study,  for  the  first  time,  what  happens  to  nanowire  conductance  as 
the  dimension  of  the  nanowire  are 
systematically  varied.  Figure  4  shows  a 
parallel  set  of  nanowire  channel  arrays, 
used  to  increase  the  sensing  area  for  greater 
sensitivity.  The  advantage  of  top-down 
lithography  based  methods  is  evident,  in 
contrast  to  bottom-up  approaches.  Other 
methods  that  do  not  use  lithography  are  not 
able  to  do  this  in  a  systematic  way.  Our 
method  has  complete  control  of  the 
geometry.  Such  studies  answer  the  question 
about  ‘why  nanotechnology’.  We  have 
shown  that  as  the  nanowire  width  is  made 
smaller  than  about  100  nm,  there  is 
dramatic  sensitivity  to  the  voltage  on  the 
top  gate.  Binding  of  the  surface  markers 
changes  the  surface  potential  and  there  is 
also  a  similar  increase  in  sensitivity  to  the 
surface  binding.  Unexpectedly,  we  found 
that  there  was  greater  sensitivity  in  the 
negative  V*  regime  (drain- source  voltage). 

This  novel  phenomenon  is  being 
investigated  for  its  greater  potential 
sensitivity.  We  found  that  there  is  a 
dramatic  increase  in  the  sensitivity  to  the 
top  gate  votage  as  the  width  of  the 
nanowires  is  decreased  below  about  100 
nm.  Figure  5  shows  that  we  have  been  able 
to  fabricate  arrays  of  50  nm  wide 
nanochannel  wires  in  a  reproducible 
manner.  Finally,  Figure  6  shows  an 
assembled  nanochannel  sample  sensor.  A  plastic  cover  that  incorporates  fluidic  channels  for  automated  sample 
exchange  using  a  syringe  pump  is  shown.  The  device  size  is  about  2  cm  in  width. 


Picture  of  nanowire  before  growing  Al  2O3 
Picture  taken  by  scanning  electron  microscope  (SEM) 


Figure  4.  Configuration  of  a  gated,  functionalized  nanowire  array. 


Fig.  5  Array  of  nanochannel  sensors 
only  50  nm  in  width. 


iBGy 


Fig.  6.  Prototype  chip  device  for 
breast  cancer  biomarker  sensing 
showing  multiple  connections,  and 
a  plastic  cover  incorporating  a 
fluidic  channels  for  automated 
exchange  of  samples. 
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We  have  now  shown  the  utility  of  the  sensors  for  detecting  a  wide  range  of  ionic  species  and 
biomarkers.  As  an  example,  we  have  shown  that  the 
nanochannel  sensors  can  sense  glucose,  urea  and  selected 
enzyme  biomarkers,  as  described  in  the  publications  in  the 
references.  The  most  important  experiments  involve  the 
testing  of  the  nanowires  against  a  selected  biomarker  target. 

We  have  shown  for  the  first  time  that  silicon  nanowires  can 
detect  the  presence  of  the  cancer  antigen,  the  mucin  related 
antibody  CA15.3.  We  examined  the  breast  cancer  serum 
biomarker  protein  CA15.3,  a  product  of  the  MUCl  gene. 

Serum  levels  of  CA15.3  reflect  tumor  burden,  and  serial 
assays  of  CA15.3  can  help  assess  disease  progression  and 
treatment  response.  The  upper  limit  of  normal  levels  of 
CA15.3  is  ~  40  Units/milliliter  (U/ml).  We  demonstrate  the 
ability  of  our  nanoscale  BIOFET  device  to  detect  CA15.3 
well  below  this  level.  Shown  is  Figure  7a  are  the  measured 
changes  in  conductance  as  the  concentration  of  the  CA15.3  is 
varied,  under  different  bias  conditions.  Figure  7b  shows  the 
differential  conductance  change  induced  as  the  antibody 
functionalized  nanowire  is  challenged  by  varying 
concentrations  of  CA15.3,  along  with  a  measurement  using 
antibiotin  as  control  (Fig.  7b).  The  data  indicate  that  the 
conductance  change  observed  is  primarily  due  to  specific 
binding  of  the  biomarker.  The  inset  in  Fig.  7a  shows  the  real¬ 
time  signal  when  exposed  to  a  calibrated  CA15.3 
concentration  of  36  U/ml.  The  device  response  time  is  less 
than  60  seconds.  The  detection  limit  of  CA15.3  concentration 
is  ~2  U/ml  (-10  pM).  Furthermore,  the  BIOFFT  response  as  a 
function  of  CA15.3  concentration,  demonstrated  here  over  a 
span  of  0  U/ml — 2,300  U/ml,  underlines  the  wide  dynamic 
range,  which  can  be  applied  in  a  wide  range  of  clinical 
situations,  for  both  early-stage  diagnosis  and  later-stage 
prognosis  without  reconfiguration.  Figure  8  shows  the  level  of 
specificity.  The  signal  response  to  CA15.3  to  a  model  protein 
shows  the  specificity.  Some  residual  non-specific  binding 
does  exist,  but  at  a  level  that  can  be  distinguished  easily  from 
background.  We  have  also  tested  the  nanochannel  sensors 
against  Prostate  Specific  Antigen  with  similar  results. 

However,  the  small  but  persistent  slope  seen  with  the 
antibiotin  challenge  took  longer  to  solve  than  expected.  We 
have  traced  the  problem  to  the  surface  functionalization 
protocol,  which  still  leaves  many  surface  binding  sites  to  ionic 
species  available.  Because  of  this,  even  though  the  response  to 
the  desired  breast  cancer  biomarker  was  clearly  greater  than 
that  for  controls,  we  postponed  studies  on  clinical  specimens. 

By  the  time  the  problem  was  resolved,  by  altering  our  surface 
surface  functionalization  protocol,  our  funding  for  the  graduate  students  working  on  the  project  had  run  out  and 
we  were  not  able  to  recruit  patients  for  studies  using  our  nanosensors.  We  regard  this  aspect  as  a  failure  on  our 
part.  On  the  positive  side,  the  nanchannel  sensors  developed  here  have  proven  to  be  of  sufficient  interest  that  a 
start-up  company,  Ninth-Sense  Inc,  is  negotiating  for  the  license  to  commercialize  the  sensors  specifically  for 
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breast  cancer  studies.  We  are  optimistic  that  the  net  result  of  our  research  will  prove  to  be  beneficial  to  the 
breast  cancer  community. 


Conclusions 

To  summarize  our  key  research  accomplishments,  we  have  developed  new  nanosensors  and  have  shown  that 
the  silicon  nanowires  are  sensitive  to  the  breast  cancer  antigen  CA15.3.  We  have  discovered  a  potentially 
extremely  interesting  “Gate  Control”  method  for  enhancing  signal  and  also  for  control  of  the  local  pH  at  the 
femtoliter  level.  We  have  also  found  a  novel  biasing  regime  which  enhances  the  sensitivity  to  a  point  where  the 
presence  of  the  biomarker  at  concentrations  well  below  clinical  levels  can  be  readily  detected.  Apart  from 
numerous  publications,  we  hope  that  the  interest  shown  in  the  technology  by  a  start-up  company  will  lead  to 
benefits  for  the  breast  cancer  community. 
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The  authors  demonstrate  the  operation  of  a  nanoscale  field-effect  pW  sensor  engineered  from  a 
functionalized  silicon  nanowire.  With  this  nanofabricated  pH  sensor,  the  change  in  the  hydrogen  ion 
concentration  or  the  pH  value  of  a  solution  can  be  detected  by  the  corresponding  change  in  the 
nanowire  differential  conductance  with  a  resolution  of  ±5nS/pH.  Fabrication  of  selective  side  gates 
on  the  nanowire  sensor  allows  field-effect  control  of  the  surface  charge  on  the  nanowire  by 
controlling  the  accumulation  of  charge  carriers  with  the  side-gate  voltage.  A  simple  physical  model 
is  used  to  analyze  the  observed  data  and  to  quantify  the  dependence  of  the  conductance  on  pH.  The 
development  of  a  nanoscale  sensor  with  physically  engineered  gates  offers  the  possibility  of  highly 
parallel  labeling  and  detection  of  chemical  and  biological  molecules  with  selective  control  of 
individual  array  elements.  ©  2006  American  Institute  of  Physics.  [DOI:  10.1063/1.2392828] 


Ultrasensitive  detection  of  biological  and  chemical  spe¬ 
cies  is  fundamental  to  the  screening  and  detection  of  disease, 
discovery  and  screening  of  drugs,  as  well  as  gas  detection 
and  biomolecular  analysis.  ’  In  particular,  the  ability  to  de¬ 
tect  ions  in  liquid  solutions  with  ion-selective  nanoscale 
electronic  sensors  is  attractive  to  a  number  of  fields.  For 
instance,  in  addition  to  the  detection  of  proteins,  virus,  and 
DNA,^~®  biocompatible  silicon-based  sensors  have  been  used 
for  in  vivo  neurological  studies  such  as  the  detection  of  cal¬ 
cium  and  potassium  ion  concentrations.  Furthermore,  highly 
parallel  detection  by  an  array  of  such  sensors  with  control 
over  individual  elements  offers  the  possibility  of  advanced 
engineering  for  enhancement  in  detection  sensitivity  as  well 
as  simultaneous  analysis  of  multiple  species. 

Nanostructures  such  as  nanowires,  nanotubes, 
nanocrystals,^  nanocantilevers,**^  and  quantum  dots"  are  par¬ 
ticularly  attractive  as  biosensors,  since  the  critical  dimen¬ 
sions  of  the  nanostructures,  such  as  the  diameter  of  the  nano¬ 
wire,  are  comparable  to  the  sizes  of  biological  and  chemical 
species.  The  detection  sensitivity  is  therefore  greatly  en¬ 
hanced  as  the  signal  can  be  effectively  transduced  because  of 
large  surface-to-volume  ratio.  In  a  nanoscale  conductor,  the 
surface-to-volume  ratio  is  large  because  of  size,  so  its  elec¬ 
trical  property  such  as  conductance  is  dominated  by  surface 
contributions.  Therefore,  the  presence  of  charged  proteins  on 
the  surface  of  an  active  nanowire  induces  a  large  fractional 
change  in  the  nanowire  conductance  and  enables  relatively 
easy  detection. 

A  silicon  nanowire  can  be  used  as  an  ultrasensitive  de¬ 
tector  by  taking  advantage  of  the  field  effect.  A  conventional 
silicon  nanowire  can  be  used  as  a  sensor  of  hydrogen  ion 
concentration  or  a  pH  sensor  by  modifying  its  surface  with 
3-aminopropyltriethoxysilane  (APTES),  which  produces 
amino  groups  as  well  as  silanol  (Si-OH)  groups  on  the  nano¬ 
wire  surface.  These  groups  operate  as  receptors  of  hydrogen 
ions,  which  undergo  protonation/deprotonation  reactions.  In 
the  process,  the  surface  charge  on  the  silicon  nanowire 
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changes,  which,  in  turn,  changes  the  nanowire  conductance. 
In  an  n-type  silicon  nanowire,  with  increasing  pH,  there  is  an 
increase  in  the  negative  charge  on  the  surface,  which  acts 
like  a  negatively  charged  gate.  This  causes  the  channel  of  the 
charge  carrier  to  deplete,  therefore  the  net  conductance  de¬ 
creases.  The  accumulation  of  carriers  by  the  reception  of 
hydrogen  ions  can  be  viewed  as  a  field  effect  as  it  modulates 
the  n-type  field-effect  transistor  (FET)  device. 

In  a  pioneering  experiment,  Cui  et  al.  have  demon¬ 
strated  the  operation  of  a  nanoelectronic  pH  sensor  fabri¬ 
cated  with  a  “bottom-up”  method.  The  change  in  the  nano¬ 
wire  conductance  as  a  function  of  pH  was  found  to  be 
consistent  with  the  previous  measurements  of  pH-dependent 
surface  charge  density  on  silica.  The  silicon  nanowires 
could  be  further  aligned  by  a  vapor-liquid-solid  growth  flow 
technique  along  the  electrodes.  Further  experiments  dem¬ 
onstrated  the  ability  of  such  sensors  to  detect  biological  spe¬ 
cies  such  as  proteins^  and  virus  molecules.^  Additional  con¬ 
trol  on  the  conductance  of  the  functionalized  FET  devices 
was  obtained  with  the  application  of  a  gate  voltage  to  the 
substrate. 

In  a  fundamentally  different  approach  to  device  architec¬ 
ture,  here  we  demonstrate  the  operation  of  a  nanoelectronic 
pH  sensor  fabricated  with  a  “top-down”  method.  The  essen¬ 
tial  advantage  of  this  approach  is  the  complete  control  over 
physical  and  electronic  degrees  of  freedom.  The  geometry 
and  alignment  of  the  nanowire  can  be  fully  controlled  by 
e-beam  lithography  and  standard  semiconductor  processing 
techniques.  Eurthermore,  physical  gate  electrodes  next  to  the 
nanowire  can  be  fabricated  with  complete  control  over  their 
location  and  size.  These  local  gate  electrodes  enable  con¬ 
trolled  accumulation  or  depletion  of  surface  charge  carriers 
on  the  nanowire  and  provide  the  ability  to  tune  the  nanowire 
conductance  necessary  for  the  optimization  of  the  detection 
sensitivity.  In  this  configuration,  the  nanowire  pH  sensor  re¬ 
gains  the  control  and  the  benefits,  usual  in  standard  elec¬ 
tronic  FET  devices. 

The  engineering  of  our  device  consists  of  two  fundamen¬ 
tal  steps;  fabrication  and  functionalization.  The  silicon  nano- 
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FIG.  1 .  (Color  online)  Device  schematic  diagram,  scanning  electron  micro¬ 
graphs.  (a)  The  schematic  diagram  of  the  silicon  nanowire  with  side  gates 
and  electrodes.  The  nanowire  is  exposed  on  three  sides  along  the  longitudi¬ 
nal  direction,  (b)  The  nanowire  shown  here  is  300  nm  wide,  230  nm  thick, 
and  8  fim  long,  (c)  A  silicon  nanowire  with  a  Au/Ti  side  gate,  (d)  The 
scanning  electron  micrograph  displays  three  silicon  nanowire  devices  on  the 
same  chip,  (e)  An  optical  micrograph  shows  the  flow  chamber  sealed  on  top 
of  the  devices  on  the  interface  board. 

wires  along  with  the  side  gates  and  the  electrodes  are  fabri¬ 
cated  by  standard  electron-beam  lithography  and  surface 
nanomachining.  The  starting  silicon-on-insulator  (SOI)  wafer 
has  a  device-layer  thickness  of  230  nm  and  oxide  layer 
thickness  of  370  nm  with  a  starting  device-layer  volume  re¬ 
sistivity  of  10-20  fl  cm.  The  device-layer  resistivity  can  be 
further  controlled  by  doping  the  wafer  by  ion  implantation  of 
boron  with  a  concentration  of  1  X  lO'^/cm^.  After  patterning 
the  nanowires  and  the  electrodes  in  separate  steps  with  sepa¬ 
rate  masks,  the  structure  is  etched  out  with  an  anisotropic 
reactive-ion  etch.  This  process  exposes  the  three  surfaces  of 
the  silicon  nanowire  along  the  longitudinal  direction  as 
shown  in  the  schematic  diagram  in  Fig.  1(a).  After  the  fab¬ 
rication  of  the  silicon  nanowire  and  the  gold  electrodes  and 
gates,  a  protective  layer  of  polymethylmethacrylate  (PMMA) 
is  spun  on  the  surface  and  only  the  silicon  nanowire  is  ex¬ 
posed  by  a  secondary  e-beam  exposure,  while  the  device 
floor  of  oxide  remains  covered.  This  process  allows  exposure 
of  only  the  silicon  nanowire  to  air/solution.  The  functional¬ 
ization  of  the  nanowire  surface  is  done  by  the  application  of 
a  2%  APTES  solution  of  methanol  for  3  h.  After  multiple 
rinsing  of  the  device  by  methanol,  the  device  is  dried  by 
nitrogen  gas  and  baked  at  80  °C  in  an  oven  for  10  min. 
Following  this  APTES  functionalization  technique,  the  plas¬ 
tic  flow  chamber  is  attached  to  the  device  with  the  applica¬ 
tion  of  PMMA  and  silicone  gel  [Pig.  1(e)].  The  flow  chamber 
is  designed  to  include  inlet  and  outlet  tubes,  connected  to  a 
syringe  pump,  to  allow  solution  flow  over  the  functionalized 
nanowire  inside  the  chamber.  The  volume  of  the  flow  cham¬ 
ber  is  estimated  to  be  —35  /xl. 

Measurement  of  pH-dependent  conductance  in  presence 

of  solution  is  carried  out  in  two  separate  methods.  One 
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FIG.  2.  Zero-bias  differential  conductance  {dl/dV)  measurements  in  solu¬ 
tion  with  different  pH  values,  (a)  dl/dV  versus  time  for  sample  1,  (b)  dUdV 
versus  pH.  (c)  Zero-bias  differential  conductance  (dl/dV)  with  increasing 
pH  values  of  the  solution  at  different  gate  voltages  for  sample  2.  (d)  Effec¬ 
tive  surface  potential  change  at  different  gate  voltages.  The  data  points  for 
Vg=0.6  V  are  shifted  by  40  nS  for  clarity  in  Figs.  2(c).  The  measurements 
are  done  with  an  ac  drive  at  37  Hz  with  an  amplitude  of  50  mV. 

method  is  the  standard  two-probe  I-V  measurement  of  the 
nanowire  along  with  a  gate  voltage  applied  to  one  of  the  side 
gates  shown  in  Fig.  1(c).  Because  of  the  high  resistance  of 
the  nanowire,  the  I-V  characteristic  is  measured  with  a  Kei- 
thley  2400  source  meter  with  a  current  resolution  of  10  pA. 
Additionally,  an  Agilent  4339B  high-impedance  bridge  is 
used  for  calibration  and  comparison.  In  order  to  discern 
small  changes  in  the  pH-dependent  conductance,  we  have 
used  another  more  elaborate  technique  of  differential  con¬ 
ductance  (dll  dV)  measurement.  The  measurement  circuit  in¬ 
cludes  a  small  ac  modulation  (provided  by  an  EG&G  5210 
lock-in  amplifier),  superimposed  on  the  dc  bias  across  the 
nanowire  (provided  by  the  Keithley  2400  source  meter) .  The 
ac  modulation  and  the  dc  bias  are  added  by  a  noninverting 
summing  circuit,  which  is  integrated  with  the  preamplifier 
circuit.  The  circuit  is  then  put  in  a  rf-shielded  aluminum  box 
to  prevent  noise  pickup.  Differential  conductance  measure¬ 
ments  are  done  by  sweeping  the  dc  bias  at  a  constant  ac 
modulation  amplitude  and  measuring  the  response  with  the 
lock-in  amplifier,  referenced  to  the  ac  signal  frequency.  This 
approach  allows  the  measurement  of  zero-bias  conductance 
dI/dV\y^Q,  avoiding  bias-dependent  effects  in  the  solution 
such  as  electrolysis. 

Figure  2  displays  zero-bias  differential  conductance  of 
the  functionalized  nanowire  in  pH  solution.  The  wire  is 
300  nm  wide,  230  nm  thick,  and  8  fim.  Solutions  with  vary¬ 
ing  hydrogen  ion  concentration  or  pH  are  made  from  phos¬ 
phate  buffered  saline,  containing  10  mM  of  phosphate  and 
130  mM  of  NaCl.  The  zero-bias  conductance  of  the  nano¬ 
wire  increases  with  increasing  pH.  Figures  2(a)  and  2(b) 
show  an  almost  linear  relation  between  dll dV  and  pH  value 
of  the  solution.  Further  evidence  of  field  effect  is  demon¬ 
strated  by  zero-bias  conductance  measurement  at  different 
gate  voltages,  shown  in  Fig.  2(c).  All  the  measurements  are 
done  at  room  temperature.  The  zero-bias  conductance  is 
tuned  by  the  application  of  the  side-gate  voltage.  A  positive 

gate  bias  implies  opening  of  the  charge  carrier  channel,  simi- 
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lar  to  the  inversion  layer  in  silicon  devices.  A  negative  gate 
bias  implies  depletion  or  squeezing  of  the  charge  carrier 
channel,  similar  to  depletion  layer.  The  zero-bias  conduc¬ 
tance  increases  or  decreases  depending  on  positive  or  nega¬ 
tive  gate  bias,  respectively. 

Field-effect  transistors  constitute  a  rich  area  of  semicon¬ 
ductor  physics.  The  term  “FET”  covers  a  large  family  of 
devices,  which  include  metal-oxide  FET  (also  called  insu¬ 
lated  gate  field-effect  transistor),  junction  field-effect  transis¬ 
tors  (JFETs),  and  so  on.  It  is  not  clear  which  model  best 
describes  the  operation  of  gated  silicon  nanowires  in  solu¬ 
tion.  The  detailed  description,  as  well  as  the  functional  de¬ 
pendence  of  device  characteristics  on  external  controls,  is 
different.  It  is  advantageous  to  describe  the  operation  of  the 
nanowire  device  in  terms  of  directly  measurable  experimen¬ 
tal  quantities  of  interest. 

Our  model  for  interpreting  the  nanoscale  devices  in  this 
study  is  based  on  a  JEET,  but  with  the  device  characteristics 
that  depend  on  solution  pH  and  ionic  strength.  Larger  scale 
devices  used  for  ion  sensitive  measurements  are  called  “ion- 
selective  field-effect  transistors”  (ISFETs).^"^  An  explicit 
comparison  with  ISFETs  may  be  made  by  considering  the 
threshold  voltage,  which  determines  the  current  flow  in  the 
silicon  wire  as^"* 

V,=  Va-^,  (1) 

where  T''  is  the  effective  surface  potential  which  is  a  function 
of  pH,  the  constant  the  ionic  charge  q 

A'F  =  -2.3«— ApH^k.  (2) 

q 

The  sensitivity  factor  a  depends  on  the  geometry  and  nature 
of  the  device  and  environmental  factors  such  as  temperature. 
It  ranges  between  0  and  I.  Equation  (2)  suggests  that  the 
temperature  dependence  of  the  surface  potential  is  compli¬ 
cated.  Eurthermore,  large  changes  in  temperature  also 
strongly  influence  transport  properties  of  the  semiconductor. 
A  detailed  study  of  temperature  effects  needs  to  take  into 
account  the  effect  on  the  concentration  and  mobility  of 
charge  carriers  in  the  semiconductor,  the  effect  of  phonons. 
Eor  ion  sensing  in  aqueous  solvents,  we  restrict  our  study 
here  to  room  temperature  (7’=  298  K).  Under  these  condi¬ 
tions,  when  a'=l, 

A^P  kT 

- - =  -  2.3a—  =  ~  -  60  mV.  (3) 

q 

For  our  device,  we  consider  the  conduction  channel  to  be 
the  central  region  of  the  nanowire  with  a  depletion  region 
between  the  channel  and  the  naturally  formed  Si02  surface. 
When  pH  of  the  solution  is  changed,  protonation  or  depro¬ 
tonation  of  the  APTES  modified  surface  can  change  the 
thickness  of  the  depletion  region.  So  the  conductance  of  the 
wire  changes.  This  is  similar  to  the  operation  principle  of 
JEET.  However,  the  wire  is  in  nanoscale.  Higher  surface-to- 
volume  ratio  of  the  smaller  wire  can  give  better  sensitivity. 
Adding  different  gate  voltages  shares  the  same  principle  of 
extending  or  squeezing  the  conduction  channel.  So  the  sen¬ 
sitivity  can  be  improved  by  adding  suitable  gate  voltage. 


In  the  differential  conductance  measurement  at  zero  bias, 
the  conductance  change  due  to  pH  is  5  nS/pH  at  Vg=0  V, 
and  the  conductance  change  due  to  gate  voltage  is  130  nS/V 
at  pH  =  8.  After  calibration  from  these  quantities,  effective 
surface  potential  change  A"*?  is  20mV/pH  at  yj=-0.4  V 
and  60  mV/pH  at  F^=0.6  V  as  shown  in  Eig.  2(d). 

Further  enhancement  of  sensitivity  of  our  sensors  can  be 
done  by  reducing  the  dimensions  of  the  silicon  nanowire  to 
increase  the  effective  surface-to-volume  ratio,  reducing  the 
doping  concentration  of  the  starting  SOI  wafer  for  optimized 
conductance,  creating  a  better  Ohmic  contact  between  the 
electrodes  and  the  underlying  silicon  surface  by  local  doping, 
and  optimizing  the  flow-chamber  design  for  faster  through¬ 
put. 

In  conclusion,  we  demonstrate  fabrication,  functionaliza¬ 
tion,  and  operation  of  a  nanoelectronic  field-effect  pH  sensor. 
The  physically  engineered  silicon  nanowire  with  side  gates  is 
fabricated  with  standard  semiconductor  processing  tech¬ 
niques.  The  functionalized  silicon  nanowire  can  be  con¬ 
trolled  with  local  nanoscale  side  gates  to  induce  inversion  or 
depletion  layers.  Our  approach  offers  the  possibility  of 
highly  parallel  detection  of  ion  or  charged  protein  and  DNA 
with  local  control  of  individual  elements.  By  selective  gat¬ 
ing,  individual  nanowires  in  an  array  can  be  turned  on  or  off 
during  functionalization.  Therefore  the  array  can  contain 
multiple  receptors  for  the  simultaneous  detection  of  multiple 
chemical  and  biological  species  in  a  single  integrated  chip. 
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Dynamical  response  of  nanomechanical  cantilever  structures  immersed  in  a  viscous  fluid  is  important  to 
in  vitro  single-molecule  force  spectroscopy,  biomolecular  recognition  of  disease-specific  proteins,  and  the 
study  of  microscopic  protein  dynamics.  Here  we  study  the  stochastic  response  of  hiofunctionalized 
nanomechanical  cantilever  beams  in  a  viscous  fluid.  Using  the  fluctuation-dissipation  theorem  we  derive 
an  exact  expression  for  the  spectral  density  of  displacement  and  a  linear  approximation  for  resonance 
frequency  shift.  We  find  that  in  a  viscous  solution  the  frequency  shift  of  the  nanoscale  cantilever  is 
determined  by  surface  stress  generated  by  biomolecular  interaction  with  negligible  contributions  from 
mass  loading  due  to  the  biomolecules. 


DOI:  10.1 103/PhysRevLett.96. 186 105 

Prom  single -molecule  force  spectroscopy  [1]  to  biomo¬ 
lecular  recognition  of  disease-specific  proteins  such  as 
cancer  antigens  [2],  micron-sized  cantilevers  have  proved 
to  be  fundamental  to  the  ultrasensitive  detection  of  small 
forces.  Usually,  forces  are  detected  by  measuring  the  de¬ 
flection  of  the  cantilever.  In  the  dynamic  case,  shift  in  the 
resonance  frequency  of  the  cantilever  is  used  to  infer  the 
magnitude  of  force.  Micromachining  techniques  now  en¬ 
able  commercial  production  of  such  cantilevers  with  di¬ 
mensions  on  the  order  of  100  /rm  as  well  as  their  routine 
use  in  force  spectroscopy. 

Decreasing  the  cantilever  dimensions  to  submicron  or 
nanometer  scales  increases  resonance  frequency  to  the 
megahertz-gigahertz  range.  The  resultant  increase  in  the 
dynamic  range  and  the  measurement  speed  can  provide  a 
better  tool  for  probing  single  molecules.  This  could  be  also 
used  for  more  sensitive  bioimaging  techniques  and  moni¬ 
toring  real-time  binding  kinetics  of  ligand-protein  binding. 
Por  biomolecular  recognition  in  a  viscous  fluid,  force  sen¬ 
sitivity  can  be  increased  by  decreasing  the  effective  viscous 
damping.  Nanoscale  cantilevers  are  hence  expected  to  have 
dramatically  enhanced  force  sensitivity  as  smaller  canti¬ 
levers  have  lower  viscous  damping. 

In  spite  of  the  importance  of  nanomechanical  cantilevers 
for  ultrasensitive  in  vitro  force  detection,  there  is  no  widely 
accepted  description  that  relates  resonance  frequency 
change  to  concentration  or  mass  loading,  over  the  entire 
range  of  viscosity  relevant  to  biomolecular  recognition  in 
viscous  fluids.  In  this  Letter,  starting  with  a  model  pro¬ 
posed  by  Sader  [3],  we  describe  the  hydrodynamics  of  a 
beam,  immersed  in  a  viscous  fluid,  while  one  of  its  sides  is 
entirely  biofunctionalized.  We  generalize  this  model  by 
including  an  additional  term  that  takes  into  account  the 
surface  stress  induced  by  the  layer  of  mutually  interacting 
biomolecules  trapped  on  the  functionalized  part  of  the 
beam  and  derive  its  exact  thermal  vibrations.  We  And 
that,  in  air,  frequency  shift  of  the  first  nanoscale  cantilever 


PACS  numbers:  81.07.-b,  45.10.-b,  82.70.Dd,  83.10.Mj 

mode  is  primarily  determined  by  surface  stress,  while 
mass-loading  effects  become  relevant  for  higher  order 
modes.  More  importantly,  in  a  viscous  solution  such  as 
water,  frequency  shift  is  dominated  by  surface  stress  rather 
than  biomolecular  mass  loading,  contrary  to  conventional 
expectation. 

Model.  — Neglecting  rotatory  inertia,  shear  deformation, 
and  internal  damping,  the  equation  of  motion  for  the  de¬ 
flection  y{x,  t)  of  a  beam  with  length  L,  width  b,  and 
thickness  d,  immersed  in  a  fluid  at  temperature  T  and 
loaded  by  a  constant  axial  force  S,  is  given  by  [3,4] 

d^y  d^v  d^v 

^  M U) =  fhix,  t)  +  fihix,  t).  (1) 

dx  dx  dr 

E  and  I  are  the  Young’s  modulus  and  moment  of  inertia  of 
the  (coated)  beam,  respectively.  Thermal  vibrations  of  the 
beam  are  induced  by  the  Brownian  (Langevin)  force  per 
unit  length  /th(x,  t).  They  are  small  enough  for  nonlinear 
convective  inertial  effects  in  the  fluid  to  be  neglected  and 
for  the  hydrodynamic  loading  /h(x,  t)  to  be  linear  in  the 
beam  displacement  [3].  The  linear  mass  (mass  per  unit 
length)  of  the  system  /x(x)  consists  of  the  linear  mass  of  the 
beam  /U^,  and  the  linear  mass  of  the  trapped  biomolecules 
pLiix).  The  axial  load  S  introduced  in  Eq.  (1)  describes  the 
mutual  interaction  of  biomolecules  adsorbed  on  the  beam 
[5].  The  boundary  conditions  for  Eq.  (1)  are  given  by 
y(0,  t)  =  y'iO,  t)  =  0,  y"(L,  t)  =  0,  and  EIy'"iL,  t)  = 
Sy'{L,  t),  where  primes  denote  spatial  derivatives  [4]. 

At  higher  concentrations,  the  biomolecules  form  a  uni¬ 
form  layer  with  linear  mass  /x/  and  thickness  h  so  that 

/jb{x)  =  /bij,  +  /bii  =  constant.  (2) 

The  mutual  interaction  of  biomolecules  within  the  layer  is 
modeled  by  taking  into  account  the  stress  a  they  generate 
on  the  coated  surface  of  the  beam.  As  shown  in  Ref.  [2], 
this  stress  enables  the  bending  of  the  silicon-nitride  micro¬ 
cantilevers  with  length  to  thickness  ratio  E/d,  ranging 
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from  10^  to  10^.  The  resulting  static  deflection,  on  the 
order  of  a  few  tenths  of  microns,  is  related  to  surface  stress 
hy  Stoney’s  formula  [6].  However,  for  the  silicon  nano¬ 
mechanical  cantilevers  under  investigation  here  {Ljd  — 
50),  Stoney’s  formula  typically  yields  angstrom-level 
hendings  (10“^  E).  So  these  nanocantilevers  remain  al¬ 
most  straight  under  the  influence  of  surface  sfress. 
Nevertheless,  as  shown  in  [5],  fhis  sfress  induces  an  effec- 
five  axial  load,  S  =  crL,  fhaf  musf  he  included  in  Eq.  (1).  In 
vacuum,  such  a  model  has  been  sfudied  in  Ref.  [7]. 

In  addifion,  hiomolecular  inferacfion  on  fhe  surface 
resulfs  in  an  effeclive  Young’s  modulus  of  fhe  layer  E/ 
[see  Ref.  [8]  ]: 

El  =  Ehh  +  Eili  ^  E^bcfilll  +  EihhSjA.  (3) 

Efjlf,  and  E//;  are  fhe  respecfive  bending  rigidifies  of  fhe 
beam  and  fhe  layer  (fhe  lasf  equably  holds  when  h  d). 
In  fhe  opposite  limil,  biomolecules  wilh  mass  m  sparsely 
scattered  over  fhe  beam  al  localions  x,  resull  in 

/r(x)  =  jjb I,  +  fiiix)  =  /uii,  +  tn^8{x  —  x,),  x,-  G  [0,  E]. 

i 

(4) 

If  fhe  average  spacing  is  large  compared  lo  Iheir  size,  Iheir 
mulual  inferacfion  is  negligible  and  5  =  0.  Considering 
lhal  Iheir  presence  does  nol  subslanlially  affecl  fhe  momenl 
Ih  of  fhe  beam,  fhe  bending  rigidily  of  fhe  whole  system  is 
fhe  same  as  for  an  unloaded  beam,  El  =  Ehlj,. 

Equations  of  Motion. — To  solve  Eq.  (1),  we  expand  fhe 
defieclion  y(x,  t)  and  fhe  force  densilies  /h(x,  t)  and  fi^{x,  t) 
in  terms  of  fhe  modes  of  fhe  bare  beam,  defined  as  fhe 
beam  wilhoul  fhe  added  mass  [/r;(x)  =  0]  Ihough  if  in¬ 
cludes  fhe  tension  5  =  aL: 

00 

y(x,  0  =  X  i  =  xlL.  (5) 

n=\ 

Similar  expressions  hold  for  fhe  force  densilies.  The  eigen- 
modes  salisfy  fhe  following  condilions: 

(p'f'if)  -  a  =  aL'^/EE, 

(6) 

4>„{0)  =  4>'„{Q)  =  0"(1)  =  0,  =  a4>'„{l). 

The  self-adjoinlness  of  Eq.  (6)  makes  fhe  modes  orthonor¬ 
mal:  pQ  4>n{^)4>i{^)df  =  8„  i.  Erom  Eq.  (6),  fhe  eigenval¬ 
ues  are  fhe  successive  posilive  roofs  of 

1  +  (1  +  £^)  coshA^j'  cosA“  -f  £„  sinhA^  sinA“  =  0,  (7) 

where  £„  =  and  Xf  =  +  si  ±  £„]'/2. 

Note  lhal,  for  £„  =  af{2p^)  1,  Eq.  (7)  reduces  lo  fhe 

usual  clamped-free  equalion,  1  -I-  cos/S„  cosh/S„  =  0.  As 
/3„  oc  „  for  large  n,  eigenvalues  for  which  n  »  ^/a  are 
essenlially  independenl  of  fhe  surface  sfress. 

Eel  /r(x)  =  /uii,  +  /xiix).  Using  Eq.  (5),  Eq.  (1)  reduces  lo 
00 

My„it)  +  k„y„it)  +  ^  0„^y/r)  =  E„_h(f)  +  E„^^it).  (8) 
;=i 


M  =  fjbi^L  is  fhe  mass  of  fhe  beam  and  E„  j,(th)(0  = 
E/„  h(th)(t)-  The  effeclive  sliffness  of  mode  n  is 

=  Eipt/L\  (9) 

The  real  and  symmelric  malrix  has  componenls 

=  E (10) 

If  /U;(x)  =  /jLi,  Ihen  O  =  M/1,  where  1  is  fhe  idenlily 
malrix  and  Mi  =  Liii  is  Ihe  layer  mass.  Equalion  (8) 
decouples  and  Ihe  mass  of  modes  y„  becomes  Ihe  lolal 
mass  of  Ihe  system,  M  +  Mi.  Bui  nonuniform  mass  dislri- 
bulions  as  in  Eq.  (4)  couple  Ihe  bare  modes  of  Ihe  beam. 

Taking  Ihe  Eourier  Iransform  of  Eq.  (8)  and  using  Ihe 
expression  for  Ihe  hydrodynamic  force  [9], 

Fn,hM  =  Mf(o^Ti(o)y„iw),  (11) 

where  Mf  =  jLpfb^  is  Ihe  mass  of  Ihe  fluid  loading  Ihe 
beam  and  r((u)  =  rr(w)  -I-  ir,(w)  is  a  complex  “hydro- 
dynamic  funclion’’  discussed  in  delail  in  [3],  we  oblain 

A(w)IKru))  =  |E,h(a;)).  (12) 

Kels  |u)  are  column  vectors  wilh  componenls  u,,  i  G  N. 
The  non-Hermitian  malrix  A((u)  is  given  by 

A(<u)  =  Ao(tu)  —  (13) 

where 

Ao(w)„y  =  -  o)\MfTi(o)  +  M„]}§„,- 

^nj  =  M„  =  M  +  0„„. 

Spectral  densities.  — As  Ihe  dissipative  (imaginary)  pari 
of  Ihe  hydrodynamic  funclion  is  frequency  dependenl,  we 
apply  Ihe  generalized  fiuclualion-dissipalion  Iheorem  [10] 
to  derive  Ihe  power  speclrum  malrix  of  Ihe  slochaslic 

forces  E„ti„  Spico)  =  |Eth(tu)){Eth(tu)r  (Ihe  overline  de¬ 
notes  Ihermal  averaging,  Ihe  superscripl  i  refers  to  Ihe 
speclral  density,  and  (E(hl  is  Ihe  Hermilian  conjugate  of 

I  Ah)): 

kT 

Spico)  =  — [A+(tu)  -  A(tu)]  =  2kTMfO)Tiicj)l,  (14) 
id) 

where  k  is  Ihe  Bollzmann  conslanl  and  T  Ihe  lemperalure. 
In  componenls,  Ihis  yields 

Fn.(.{co)F*p,i,{ay')  =  2kTMpcoT i{(o)8,^8{(o  -  co').  (15) 

Nolice  lhal  Ihis  expression  does  nol  depend  on  O.  Il  is  Ihe 
same  as  for  a  bare  beam.  As  seen  above,  Ihe  slochaslic 
forces  acling  on  dislincl  modes  are  uncorrelaled. 
Neverlheless,  Iheir  power  speclrum  is  nol  conslanl,  con- 
Irary  to  Ihe  assumplion  made  in  Ref.  [3],  for  Ihe  frequency 
dependence  of  Ihe  dissipalive  part  of  Ihe  hydrodynamic 
funclion  makes  Ihem  non-Markovian.  Equalion  (15)  is  Ihe 
generalizalion  of  Ihe  expression  derived  by  Paul  and  Cross 
for  a  single  canlilever  mode  [11].  Now,  inverling  Eq.  (12), 
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we  obtain 

Iy(^^))  =  Y(")l^th(«^)).  (16) 

and  |})(«))())(«)|  =  ;t'(<n)|Fth(«))(fth(«^)lY^(")-  Spectral 
averaging  the  latter  and  using  Eq.  (14),  we  find  the  power 
spectrum  matrix  of  the  deflection  modes 

5,(«)  =  2kTMfCoYi{w)x{(o)xH(^).  (17) 

Introducing  \(p^)  with  components  (f>ni^),  the  Eourier 
transform  of  the  deflection  (5)  reads  co)  =  {(t>^\y{o})) 
and  using  Eq.  (17),  we  find  its  spectral  density  to  be 


\y{x,  <u)P“  =  2kTMfa)riia)){(f)^\xio))x^io))\(f>{)-  (18) 

The  total  mass  of  the  particles  trapped  on  the  beam  Mi  is 
small  compared  to  the  mass  of  the  beam.  This  justifies 
treating  perturb atively  provided  its  elements  stay 

small  compared  to  the  diagonal  elements  of  Ao(&>).  This 
is  indeed  the  case  [12]  provided  MfTj{cj)  »  4M;.  By 
inverting  Eq.  (13),  we  obtain  in  first  order  in  O 

=  YoYo  +  +  Yo‘I>)Yo  +  0(<I>^),  (19) 

where  xo  =  *•  Reinstating  in  Eq.  (18),  we  finally  get 


\yix,  (w)P’  =  2kTM f(oT i{o)) 


L,4i«;K 


20,=  ££ 


n=\  p=l 


+  0(<I)2) 


(20) 


In  this  expression,  ^  =  xjL,  is  given  by  Eq.  (14),  and 
the  quantity  A„  =  Aq{co)„„/M„  reads 

A„  =  (ol-  co\l  +  A„r(u;)],  (21) 

where  tu„  =  is  the  frequency  in  vacuum,  A„  = 

Mf/M„,  and  9f(A„)  is  the  real  part  of  A„.  An  expression 
similar  to  Eq.  (20)  can  be  derived  for  the  slope  of  the 
deflection  provided  4>j{^)  is  replaced  by  the 

overall  prefactor  is  divided  by  L^. 

Expression  (20)  is  valid  for  any  mass  distribution  /x;(x) 
along  the  beam.  Eor  a  uniform  layer  with  linear  mass  /r,/, 
<I>„^  =  Lfjii8„p,  and  then  =  0.  All  modes  have  the 
same  effective  mass,  M„  =  M  +  Lfii,  and  are  decoupled. 
Reinstating  in  Eq.  (20),  the  second  term  vanishes  and  we 
obtain  the  exact  spectral  density  of  a  composite  beam 
consisting  of  the  original  beam  plus  the  layer.  Eor  mole¬ 
cules  trapped  on  the  beam  at  positions  f the  mass  profile 
given  in  Eq.  (4)  leads  to  and 

Eq.  (20)  is  valid  up  to  first  order  in  m  provided  the 
frequency  satisfies  MfTiio))  »  4Mp  Interestingly,  if  we 
assume  N  molecules  to  be  randomly  scattered  along  the 
beam  in  a  uniform  way  and  average  <I>„p  accordingly,  we 
find  (($„p))  =  mXi  /  =  NmS„p.  As  the 

total  mass  of  the  trapped  molecules  is  small  compared  to 
the  mass  of  the  beam,  in  the  first  approximation,  the 
average  spectral  density  is  the  same  as  the  spectral  density 
of  their  average  mass  distribution — i.e.,  the  spectral  den¬ 
sity  of  a  uniform  layer  of  mass  Nm. 

In  Eig.  I,  we  compare  the  bare  beam  (O  =  0)  spectral 
density  of  the  deflection  slope  at  the  tip  of  a  rectangular 
silicon  nanocantilever  (E  =  160  GPa,  p  =  2.33  X 
10^  kg/m^)  to  Sader’s  result  [3]  in  air  and  water.  The 
beam  dimensions  are  (iXfiXL  =  0.2X0.2X10  /rm. 
Erom  [13],  at  r  =  27  °C,  the  viscosities  are  p^ir  =  1.86  X 
10“^  "fiwater  =  8.59  X  10  ^  kg/(m.s)  and  the  densities  are 
Pair  =  1-18,  Pwater  =  997  kg/m^.  Typical  values  for  the 
Reynolds  number  defined  as  [3]  Re  =  pfb'^co/Ari  are 
Reair  —  0.01  and  Re^ater  —  0.036  at  resonance.  Although 
different,  Sader’s  formula  can  be  shown  to  reduce  to  (20) 
provided  A|r(tu)|  1.  This  explains  why  the  results  are 


very  similar  in  air  [A|r(<w)|  0.06  at  resonance]  while 

they  start  to  differ  in  water  [A|r(<w)|  25  at  resonance]. 

Frequency  shift.  — As  stated  earlier,  when  trapped  mole¬ 
cules  form  a  uniform  layer,  the  exact  spectral  density  of  the 
beam  deflection  is  given  by 


ly(Y  u;)P^ 


2kTMfCoTiico)  - 

{M  +  Mif  |A„P 


(22) 


where  M;  is  the  mass  of  the  layer  and  where  A„  is  given  in 
Eq.  (21)  with  =  M  +  Mp  When  the  peaks  of  Eq.  (22) 
are  sharp  enough,  the  hydrodynamic  function  is  almost 
constant  in  their  vicinity  and  the  resonant  frequency  sat¬ 
isfies  the  self-consistent  equation: 

,  ,  ,  1 R  +  V4R2  3/2 

^R.n=  f((^R.n,M„)k„,  f  =  - - - ’  (^3) 

where  R  =  M„  +  MfT^,  I  =  MfTi,  and  Er,  =  rr_,(wp_„). 
Erom  the  expression  (23),  the  mass  and  stiffness  variations 
due  to  the  layer,  SM  and  induce  a  relative  frequency 
shift  between  a  bare  and  a  loaded  beam: 


^^R.n 

^R,n 


2 


Sk„ 

k„ 


a  Inf 

-f - -SM 

dM 


(24) 


f=  00/271  (MHz) 

FIG.  1.  Inset:  schematic  diagram  of  the  beam  with  a  layer  of 
molecules  on  its  functionalized  side.  Main  figure:  exact  spectral 
density  of  the  deflection  slope  of  a  silicon  nanobeam  without 
molecular  layer.  Calculations  done  in  air  and  water  (solid  lines) 
compared  to  Sader’s  results  [3]  (dashed  lines). 


186105-3 


PHYSICAL  REVIEW  LETTERS 


week  ending 
12  MAY  2006 


PRL  96,  1 86 1 05  (2006) 


a  a 

FIG.  2.  Relative  frequency  shift  (in  %)  in  air  and  water  vs  the  dimensionless  stress,  a  =  crL? lEI.  Mi  =  M/100  and  £;  =  0.  Other 
parameters  are  the  same  as  in  Fig.  1.  Solid  line:  exact  result  from  Eq.  (22).  Dotted  line:  linear  approximation  from  Eqs.  (24)  and  (25). 


Here,  dM  =  Mi  and  its  prefactor  in  Eq.  (24)  takes  into 
account  dissipative  (MyrE,)  and  fluid  mass  loading 
effects.  According  to  (9),  the  two  contributions  to  the 
stiffness  8k„  come  from  the  bending  rigidity,  Ei,Ii, 
Efjlf,  +  Eili,  and  from  the  surface  stress  through  the  eigen¬ 
value  Erom  Eqs.  (3)  and  (7),  we  And 

^  3hEi  2T„t„  +  (3o„iT„  +  f„) 

K  dE,  l3lSt„-T„) 

where  /Sq  „  is  the  nth  root  of  cosh(jS)  cos(/S)  -1-1=0  and 
where  T„  =  tanhjSo,,,  t„  =  tan^So.,,.  The  last  term  of  (25) 
has  been  obtained  from  (7)  in  perturbation.  It  is  valid  when 
a  and  vanishes  as  n  — ►  oo. 

Using  the  same  data  as  in  Eig.  1,  we  display  in  Eig.  2  the 
relative  frequency  shift  in  air  (left)  and  water  (right)  versus 
the  dimensionless  surface  stress  a.  The  exact  shift  is 
evaluated  from  the  spectral  density  in  Eq.  (22)  and  com¬ 
pared  to  its  linear  approximation  in  Eq.  (24).  The  layer 
mass  has  been  arbitrarily  Axed  to  1%  of  the  beam  mass 
and  El  set  to  zero,  hence  the  negative  offset  observed  in 
air  at  O'  =  0.  Eor  typical  values  of  the  surface  stress,  a  ~ 
10“^  Jm“^  [see  Wu  et  al.  in  Ref.  [2]],  lal  ^  1.  In  air, 
A|r(<Ws„)|  I,  and  f~l/M.  Then,  dhif/dM~ 
—  XjM,  cHn//5m/j„  ~  0,  and  we  recover  the  usual  fre¬ 
quency  shift  for  a  linear  oscillator  in  vacuum.  As  seen  on 
the  left  panel,  the  first  peak  is  the  most  sensitive  to  a.  The 
deviation  of  the  data  from  the  linear  result  in  Eq.  (25) 
indicates  that  the  condition  a  with  n  =  I  becomes 
violated.  This  effect  disappears  for  the  second  and  third 
peaks  that  are  less  sensitive  to  a.  In  water  (right  panel),  a 
single  broad  peak  occurs.  Equation  (23)  loses  its  accuracy 
but  the  frequency  shift  in  Eq.  (24)  derived  from  it  is  still 
acceptable.  The  contribution  of  d  In// dM  becomes  negli¬ 
gible  while  d\nf/do)if„ - dlnEi/dtu^,,  becomes  im¬ 

portant,  hence  the  increase  in  the  slope  of  the  relative 
frequency  shift  versus  a  in  water  compared  to  air. 

In  conclusion,  we  treat  the  thermal  response  of  biofunc¬ 
tionalized  nanocantilevers  with  a  generalized  fluctuation- 
dissipation  relation.  In  a  viscous  fluid  like  water,  the  reso¬ 


nance  frequency  shift  for  a  continuous  distribution  of  bio¬ 
molecules  on  the  cantilever  surface  is  found  to  be 
dominated  by  surface  stress  rather  than  the  mass  loading 
of  biomolecules. 
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We  report  amplification  of  biomolecular  recognition  signal  in  lithographically  defined  silicon 
nanochannel  devices.  The  devices  are  configured  as  field  effect  transistors  (FET)  in  the  reversed 
source-drain  bias  region.  The  measurement  of  the  differential  conductance  of  the  nanowire  channels 
in  the  FET  allows  sensitive  detection  of  changes  in  the  surface  potential  due  to  biomolecular 
binding.  Narrower  silicon  channels  demonstrate  higher  sensitivity  to  binding  due  to  increased 
surface-to-volume  ratio.  The  operation  of  the  device  in  the  negative  source-drain  region 
demonstrates  signal  amplification.  The  equivalence  between  protein  binding  and  change  in  the 
surface  potential  is  described.  ©  2007  American  Institute  of  Physics.  [DOI:  10.1063/1.2822445] 


Ultrasensitive  recognition  of  molecules  is  important  in 
basic  science,  as  well  as  in  pharmacology  and  medicine  for 
the  analysis  of  biomarkers.  ’  Nanotechnology  has  made  it 
possible  to  enhance  detection  sensitivity  by  using 
nanowires,^’"^  nanotubes,^  nanocrystals,®  nanocantilevers, 
and  quantum  dots.^  The  signal  arising  from  biomolecular 
binding  on  the  surface  can  alter  measurable  physical  proper¬ 
ties  of  the  nanoscale  device  such  as  electrical  conductance 
with  increased  sensitivity  due  to  the  large  surface-to-volume 
ratio. 

A  silicon  nanowire  can  be  used  as  the  source-drain  chan¬ 
nel  of  a  field  effect  transistor  (FET).  In  conventional  FETs, 
lithographic  methods  are  used  to  fabricate  gates  at  the  bot¬ 
tom,  the  top,  or  the  side.  In  addition,  the  nanochannel  surface 
can  be  functionalized  with  specific  receptor  or  antibody.  In  a 
fluid,  the  ligand  (or  antigen)  can  bind  to  the  receptor,  which 
results  in  a  change  in  the  surface  charge  profile  and  the  sur¬ 
face  potential.  Essentially,  this  binding  behaves  as  a  field 
effect.  The  conductance  and  the  I-V  characteristics  of  the 
nanowire  can  therefore  be  used  to  characterize  biomolecular 
binding — for  instance,  to  determine  concentration  and  bind¬ 
ing  dissociation  constant.  In  this  letter,  we  show  that  charac¬ 
teristics  of  differential  conductance  dlldV  can  be  used  for 
even  higher  sensitivity  in  the  field  effect  due  to  biomolecular 
binding.  In  particular,  dlldV  characteristics  allow  measure¬ 
ment  at  low  bias,  essential  for  avoiding  electrolysis. 

In  a  series  of  experiments,  Lieber  and  co-workers  have 
demonstrated  the  operation  of  nanowire  FET  sensors  fabri¬ 
cated  with  a  bottom-up  method,  using  chemically  grown  sili¬ 
con  nano  wires.  These  nano  wire  field-effect  sensors  show  sig¬ 
nificant  advantages  of  real-time,  label-free  and  highly 
sensitive  detection  of  a  wide  range  of  analytes,  including 
proteins,  nucleic  acids,  small  molecules,  and  viruses  in 
single-element  or  multiplexed  formats.*®  These  bottom-up 
approaches,  however,  have  limitations  arising  from  the  lack 
of  control  in  materials  and  fabrication  engineering.  In  con¬ 
trast,  top-down  methods  of  fabrication  by  e-beam  lithogra¬ 
phy  or  submicron  optical  lithography  enable  scalable  manu¬ 
facturing. 


“'Electronic  mail:  mohanty@physics.bu.edu. 


Here,  we  demonstrate  that  the  differential  conductance 
dlldV  characteristics  of  silicon  nanowire  field  effect  transis¬ 
tors  with  width  below  200  nm  show  strong  dependence  on 
gate  voltage.  We  use  the  peak  position  of  the  dll  dV  curve  to 
characterize  antibiotin  detection  with  sub-ng/ml  sensitivityby 
functionalizing  the  silicon  nanowire  surface  with  biotin. 
From  this  measurement,  we  can  extract  the  dissociation  con¬ 
stant  of  biotin- antibiotin  binding  to  be  5.2  X  10“*®  M,  in 
good  agreement  with  other  measurements. 

By  functionalizing  the  silicon  nanowire  surface  with  bi¬ 
otin,  we  show  that  the  device  can  be  used  to  characterize 
protein  binding.  Furthermore,  we  show  equivalence  between 
the  field  effect  created  by  the  biomolecular  binding  and  the 
reference  gate  voltage.  This  equivalence  allows  the  charac¬ 
terization  of  ligand-receptor  binding  by  determining  the  cali¬ 
brated  gate  voltage  dependence  rather  than  measuring  the 
more  cumbersome  concentration  dependence. 

Our  devices  are  fabricated  by  standard  e-beam  lithogra¬ 
phy  and  surface  nanomachining.  The  starting  (100)  SOI  wa¬ 
fer  has  a  device  layer  thickness  of  100  nm  and  oxide  layer 
thickness  of  380  nm,  on  a  600  ptm  boron-doped  substrate. 
The  device-layer  volume  resistivity  ranges  from 
10-20  H  cm.  After  patterning  the  nanowires  and  the  elec¬ 
trodes  in  separate  steps,  the  structure  is  etched  out  with  an 
anisotropic  reactive-ion  etch.  This  process  exposes  the  three 
surfaces  of  the  silicon  nanowire  along  the  longitudinal  direc¬ 
tion,  resulting  in  structure  with  three-dimensional  relief  and  a 
large  surface-to-volume  ratio.  Finally,  a  layer  of  AI2O3 
(5-20  nm  thick)  is  grown  by  atomic  layer  deposition.  Aplas¬ 
tic  fiow  cell  is  used  to  bathe  the  device  in  a  small  volume  of 
fiuid  (20-30  fjuL)  avoiding  the  mixing  problem  of  microfiu- 
idic  channel.** 

The  device  pictures  and  measurement  configurations  are 
shown  in  Fig.  1.  The  use  of  parallel  nano  wires  improves  the 
average  signal  while  maintaining  a  high  surface-to-volume 
ratio.  A  fabricated  top  gate  is  used  to  characterize  current- 
voltage  relation  of  the  device  in  Fig.  1(c).  Small  changes  in 
the  conductance  of  the  device  are  best  measured  by  consid¬ 
ering  the  differential  conductance  g(U(js,yg)  =  (o'/(js/6'y;js)y 
with  the  measurements  made  at  constant  top  gate  voltage  Vg. 
In  biological  sensing  measurements  [Fig.  1(d)],  an  Ag/AgCl 
reference  gate  is  used.  The  conductance  is  a  function  of  the 
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FIG.  1 .  (Color)  Device  configuration.  Scanning  electron  micrographs  of  (a) 
multiple  sensor  array  and  (b)  single  wire,  (c)  Schematic  of  device  cross- 
sectional  view  with  top  gate  geometry,  (d)  Schematic  of  the  cross-sectional 
view  of  the  device,  showing  reference  gate  geometry. 


bias  voltage  Vjjj,  the  gate  voltage  Vg,  or  equivalently  the 
surface  potential,  which  can  be  changed  either  by  changing 
Vjg  or  the  concentration  of  charged  biomolecules. 

Figure  2(a)  is  the  characteristic  of  the  device  with  a  wire 
width  of  350  nm.  Figure  2(c)  shows  the  characteristic  I-V 
curve  for  a  80  nm  device.  In  Figures  2(b)  and  2(d),  each 
curve  corresponds  to  differential  conductance  measurements 
at  fixed  values  of  Vg,  ranging  from  -1  to  -i- 3  V.  At  large 
negative  all  the  curves  start  with  a  common  nonzero 
value  of  g.  Such  a  universal  value  suggests  that  conductance 
is  dominated  by  the  Schottky  barrier  at  the  junction.*^  Figure 
2(b)  shows  differential  conductance  for  the  350  nm  wire  as  a 
function  of  which  shows  a  peak  in  differential  conduc¬ 
tance  near  0  V.  As  F(js  is  increased  above  a  threshold  voltage 
Fth.  the  slopes  of  both  /-V  characteristic  and  differential  con¬ 
ductance  increase.  For  the  350  nm  wire,  the  differential  con¬ 
ductance  peak  positionin  Fig.  2(b)  is  insensitive  to  Vg.  In 
contrast,  the  peak  position  in  Fig.  2(d)  is  strongly  dependent 
on  Vg.  The  sensitivity  of  the  conductance  peak  of  the  device 
with  narrower  width  comes  from  the  threshold  voltage 
change.  In  our  model,  is  linearly  dependent  on  the 
surface-to-volume  ratio 


(w  +  2H)L  W  +  2H 
wHL  ~  wH 

Here,  w  is  the  width  of  the  wire,  H  is  the  thickness  of  the 
device  layer,  and  L  is  the  length  of  the  wire.  The  peak  posi¬ 
tion  of  the  differential  conductance  at  zero  top  gate  voltage  is 
plotted  as  a  function  of  the  nanowire  width  in  Fig.  2(e). 
Following  the  methods  given  in  Ref.  14,  the  conductance 
peak  of  the  dH  dV  versus  F(js  can  be  calculated  in  a  semicon¬ 
ductor  model.  With  increased  surface-to-volume  ratio,  and 
hence  a  higher  threshold  voltage,  the  conductance  channel  is 
closed  at  zero  bias.  A  higher  gate  voltage,  i.e.,  a  larger  nega¬ 
tive  bias  voltage,  is  then  needed  to  open  the  narrow  wire 
channel.  Therefore,  the  conductance  peak  shifts  to  negative 
drain  voltage  [Fig.  2(e),  red  points]. 

In  the  experiment,  the  reference  gate  is  immersed  in  the 
solution,  and  a  back  gate  voltage  is  applied  to  the  substrate. 
In  an  array  of  sensors,  the  solution  reference  gate  voltage  F^g 
is  common  to  all  the  nanowires.  However,  the  source-drain 
voltage  Fds  for  each  nanowire  in  the  nanosensor  array  can  be 
controlled  individually.  In  this  more  flexible  configuration, 
the  high  voltage  required  for  the  back  gate  can  be  avoided." 
At  the  same  time,  after  channel  is  opened  gradually  by  F(js, 
the  sensitivity  of  the  device  can  be  amplified.  This  naturally 
implies  a  region  of  operation  at  negative  bias  for  high  sensi¬ 
tivity  measurements. 

The  signal  amplification  effect  in  reverse  bias  is  demon¬ 
strated  in  protein  (antibiotin)  measurements.  Silicon  nano¬ 
wires,  functionalized  with  biotinylated  bovine  serum  albu¬ 
min,  are  used  to  detect  antibiotin  in  1  mM  NaCl  and  1  mM 
phosphate  buffer  solution.  At  the  concentration  of  salt  used 
in  solution,  the  Debye  screening  length  Xyj  at  room  tempera¬ 
ture  is  about  9.6  nm.  X.^,  is  sufficiently  large  that  the  surface 
potential  is  sensitive  to  protein  binding,  but  short  enough  to 
screen  out  biomolecules  in  solution.""^  All  solutions  used  in 
the  measurements  are  dialyzed  to  maintain  constant  ionic 
strength.  The  differential  conductance  change  Ag  due  to 
100  ng/ml  antibiotin  in  the  solution  is  0.02±0.01  nS  at  F(js 


FIG.  2.  (Color)  Width  dependence  of  the  device  I-V  characteristics,  [(a)  and  (b)]  /jj  and  dl/dV  vs  source  drain  voltage  Vjs  at  different  top  gate  voltage  Vg 
(-1  -  -1-3  Vwith  1  V  step  increase)  for  350  nm  wires,  [(c)  and  (d)]  and  dl/dV  vs  Vd.  at  different  Vg  (—1-3  V  with  1  V  step  increase)  for  80  nm  wires,  (e) 
Peak  position  of  the  dlldV  peak  at  V^=0  V  for  devices  with  different  width.  Red  points  show  the  simulation  results  used  in  the  model  in  Ref.  14.  The  inset 
is  the  differential  conductance  as  the  function  of  source  drain  voltage  at  V^=0  V  for  different  wire  widths. 
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FIG.  3.  (Color)  [(a)  and  (b)]  Differential  conductance 
vs  time.  The  phosphate  buffer  is  first  injected,  followed 
by  a  100  ng/mL  of  antibiotin  in  the  same  buffer,  with 
data  acquisition  being  stopped  for  30  s  following  the 
injection.  The  bias  is  set  at  yjjs=-0.4  V  and 
-0.9  V,  respectively,  (c)  Change  in  the  differential  con¬ 
ductance  of  the  device  introduced  by  antibiotin  at  dif¬ 
ferent  Vds  (bottom  axis,  with  fixed  yj.g=0.3  V,  black 
solid  squares).  Red  solid  triangles  are  the  data  at  differ¬ 
ent  Vjg  (top  axis,  with  Vds  =  0  V).  Inset  is  the  signal 
noise  ratio  of  the  device  at  different  V^s-  (d)  Compaii- 
son  of  conductance  change  introduced  by  5  mV  of  ref¬ 
erence  gate  voltage  change  (black  dots,  left  axis)  and 
100  ng/ml  of  antibiotin  solution  (red  triangles,  right 
axis). 


=  -0.4  V  [see  Fig.  3(a)],  while  Ag  is  45  ±0.1  nS  at  1^5= 
-0.9  V  [see  Fig.  3(b)].  The  signal-to-noise  ratio  increases 
from  2  to  34  [see  Fig.  3(c)  inset].  The  signal  due  to 
100  ng/ml  antibiotin  injection  is  plotted  as  a  function  of  V^s 
and  Fjg  [see  Fig.  3(c)].  This  plot  clearly  shows  the  effect  of 
reverse  bias  amplification. 

The  change  Ag  above,  due  to  concentration  change  at 
fixed  reference  gate  voltage,  can  be  compared  to  the  change 
caused  by  varying  the  reference  gate  voltage  while  keeping 
the  concentration  fixed.  An  equivalence  between  the  surface 
potential  change  and  concentration  change  is  then  estab¬ 
lished.  Ag  introduced  by  100  ng/ml  (—660  pM)  of  antibi¬ 
otin  is  equivalent  to  a  gate  voltage  change  of  7.2±0.3  mV 
[see  Fig.  3(d)]. 

The  fundamental  advantage  of  our  label-free  device  ar¬ 
chitecture  is  the  combination  of  high  detection  sensitivity 
and  standard  semiconductor-based  fabrication  techniques, 
suitable  for  scalable  manufacturing.  We  operate  the  device  in 
the  reverse  bias  region  without  applying  high  voltage  to  the 
back  gate'^  or  to  the  reference  gate  voltage  in  the  solution. 
Our  device  configuration  allows  device-level  signal  amplifi¬ 
cation  and  the  required  degree  of  control  to  enable  large- 
scale  parallel  architecture  for  detection  of  multiple  target 
molecules.  This  combination  is  important  for  any  clinical 
application,  including  disease  screening,  diagnosis,  monitor¬ 
ing,  and  prognosis,  all  of  which  require  fast  analysis  time, 
small  sample  volume,  and  low  cost. 

In  conclusion,  we  fabricate  silicon  nanowire  channel 
FET  sensors  by  a  top-down  e-beam  lithography  method.  Dif¬ 
ferential  conductance  characteristics  of  the  devices  show  de¬ 
pendence  on  the  source-drain  channel  width  of  the  nanoscale 
sensor.  We  take  advantage  of  the  strong  dependence  in  the 
reverse-bias  region  to  demonstrate  amplification  of  biomo- 
lecular  binding  signal  by  modulating  the  top  gate  voltage. 


The  operation  of  the  device  at  low  reference  gate  voltage 
prevents  electrolysis.  By  comparing  the  signal  generated  by 
protein  binding  with  the  change  in  the  reference  gate  voltage, 
we  establish  an  equivalence  between  the  effect  arising  due  to 
the  concentration-dependent  surface  charge  density  and  the 
tunable  surface  potential  in  the  silicon  nano  wires. 
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Silicon-based  nanochannel  glucose  sensor 
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Silicon  nanochannel  biological  held  effect  transistors  have  been  developed  for  glucose  detection. 
The  device  is  nanofabricated  from  a  silicon-on-insulator  wafer  with  a  top-down  approach  and 
surface  functionalized  with  glucose  oxidase.  The  differential  conductance  of  silicon  nanowires, 
tuned  with  source-drain  bias  voltage,  is  demonstrated  to  be  sensitive  to  the  biocatalyzed  oxidation 
of  glucose.  The  glucose  biosensor  response  is  linear  in  the  0.5-8  mM  concentration  range  with 
3-5  min  response  time.  This  silicon  nanochannel-based  glucose  biosensor  technology  offers  the 
possibility  of  high  density,  high  quality  glucose  biosensor  integration  with  silicon-based  circuitry. 
©  2008  American  Institute  of  Physics.  [DOI:  10.1063/1.2832648] 


Field  effect  devices,  such  as  capacitive  electrolyte- 
insulator-semiconductor  sensor,  light-addressable  potentio- 
metric  sensor,  and  ion-sensitive  held  effect  transistor 
(ISFET)  for  glucose  detection,*”^  have  been  extensively  stud¬ 
ied  in  recent  years.  Although  these  devices  are  limited  by  the 
dependence  of  the  sensor  response  on  buffer  capacity,  ionic 
strength,  and  /jH  of  the  test  sample,  their  compatibility  with 
advanced  microfabrication  technology  may  enable  their  po¬ 
tential  commercialization.  Glucose  biosensor  is  a  particularly 
attractive  enzyme  biosensor  due  to  its  potential  widespread 
use  in  clinical  applications. 

Currently,  glucose  detection  is  mostly  limited  to  in  vitro 
test  of  blood  samples,  although  it  is  more  meaningful  to  per¬ 
form  in  vivo  test  by  implantable  sensing  devices  for  contin¬ 
ued  monitoring  of  blood  glucose  level.  To  this  end,  nanoscale 
sensors  may  be  fundamentally  valuable.  Semiconductor 
nanowires,  grown  from  bottom-up  approach,  have  been  dem¬ 
onstrated  as  good  candidates  for  ultrasensitive  biosensors  in 
many  applications."*  ®  However,  most  of  the  existing  studies 
based  on  bottom-up  approaches  face  the  limitation  of  com¬ 
plex  integration  and  scalable  large-scale  manufacturing. 
Fabrication  of  silicon  nanoscale  devices  with  top-down 
approaches  ’  using  lithography  has  the  inherent  benefit  of 
standard  semiconductor  processes  and,  hence,  better  control 
of  device  properties.  Therefore,  field  effect  devices,  using 
nanoscale  technology,  offer  the  possibility  of  high- 
performance,  low-cost  implantable  glucose  biosensors. 

Here,  we  demonstrate  silicon  nanochannels,  surface 
functionalized  with  glucose  oxidase,  as  a  field  effect  glucose 
biosensor.  The  glucose  biosensor  with  a  set  of  silicon  nano¬ 
wires  as  nanochannels,  50- 100  nm  wide,  100  nm  high,  and 
6  /um  long,  is  fabricated  from  a  silicon-on-insulator  (SOI) 
wafer.  Figure  1  shows  a  typical  scanning  electron  micro¬ 
graph  of  the  device.  The  SOI  wafer  consists  of  a  100  nm 
thick  silicon  top  layer,  a  500  ptm  silicon  substrate,  and  a 
380  nm  Si02  insulation  layer  in  between.  The  silicon  top 
layer  is  lightly  doped  with  boron  concentration,  10“'^  cm"^, 
as  the  device  layer.  The  silicon  nanowires  are  patterned  with 
electron  beam  lithography,  and  two  sidewalls  are  exposed 
after  reactive  ion  etching.  The  Ti/Au  layer  is  deposited  with 
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thermal  evaporator  to  function  as  the  source  and  drain  con¬ 
tact  electrodes  without  further  annealing.  The  silicon  nano¬ 
wires  are  covered  with  a  layer  of  10  nm  AI2O3,  grown  by 
atomic  layer  deposition,  to  prevent  current  leakage  between 
analyte  solution  and  silicon  nanowires. 

Before  modification,  the  AI2O3  surface  is  treated  with 
oxygen  plasma^  for  two  purposes:  to  clean  the  sample  sur¬ 
faces  and  to  generate  a  hydrophilic  surface.  The  wires  are 
first  modified  by  APTES  (3%  in  ethanol  with  5%  water). 
Then  3%  glucose  oxidase  in  acetic  chloride  (50  mM)  buffer 
(/jH  of  5.1)  (5%  glycerol,  5%  BSA)  is  deposited  on  the 
sample  and  kept  in  glutaraldehyde  vapor  for  40  min.  The 
sample  is  dried  in  air  for  15  min.  Glucose  samples  are  made 
in  solution  with  50  mM  NaCl  and  50  mM  (or  81.8  mM)  of 
potassium  ferricyanide. 

All  electrical  measurements  are  done  at  room  tempera¬ 
ture.  The  voltage  V  across  the  nanowires  is  applied  by  using 
a  modulated  sine  wave  source  upon  dc  bias.  The  current 
I  across  the  nanowires  is  converted  into  voltage  by  a  feed¬ 
back  resistor,  and  this  voltage  is  detected  with  a  lock-in  am¬ 
plifier.  A  voltage  is  applied  on  the  Ag/  AgCl  reference  elec¬ 
trode  contacted  with  solution  as  the  reference  gate  voltage, 

(Ref.  7).  This  allows  direct,  real-time  measurement  of 
the  differential  conductance  g  =  {dll dV)y  .  In  particular,  the 
method  allows  studies  at  zero  bias  and  reverse  bias  as  well. 


FIG.  1.  (Color  online)  (a)  Scanning  electron  micrograph  of  silicon  nano¬ 
wires,  the  working  part  of  the  biosensor,  (b)  Top  view  of  a  50  nm  wide 
nanowire  shows  the  controlled  lineal'  geometi'y. 
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FIG.  2.  (Color  online)  (a)  The  glucose  biosensor  response  in  the  0.5-8  mM 
concentration  range,  which  is  fitted  into  a  lineai*  curve.  y(js=-0.69  V, 
\/g=0  V.  All  solutions  contain  50  mM  of  NaCl  and  50  mM  of  potassium 
ferricyanide.  The  top  inset  shows  the  real  time  differential  conductance 
change  due  to  exchange  glucose  solutions.  The  bottom  inset  is  differential 
conductance  of  the  device  versus  glucose  concentration.  The  random  en'or  is 
estimated  to  be  5%  of  the  device  response  due  to  the  addition  of  different 
glucose  solutions,  (b)  Device  responses  for  81.8  mM  (red  curve)  and 
50  mM  (black  curve)  potassium  ferricyanide  solutions  are  displayed. 

Figure  2(a)  gives  the  real-time  device  response  when 
we  add  glucose  solutions  upon  initial  use.  We  find  that  there 
is  a  linear  relation  between  device  response  and  glucose  con¬ 
centration  in  the  range  of  0.5-8  mM.  From  the  response 
curve  [the  bottom  inset  of  Fig.  2(a)],  we  find  that  the  differ¬ 
ential  conductance  change  of  the  device  reaches  saturation 
after  the  glucose  concentration  of  16  mM.  We  can  also  esti¬ 
mate  the  diffusion-limited  time  response  of  the  device  to  be 
3-5  min,  which  is  the  same  as  the  typical  response  times  of 
thick-layer  glucose  oxidase  modified  ISFET.*"  The  response 
time  could  be  further  improved  by  enzyme  monolayer 
functionalization.  *  * 

The  oxidation  of  glucose  catalyzed  by  glucose  oxidase  is 
a  two-substrate  reaction.*^  One  molecule  of  glucose  oxidase 
contains  two  molecules  of  oxidized  form  of  flavine  adenine 
dinucleotide  (FAD).  In  the  first  substrate  reaction,  there  is  a 
formation  of  an  enzyme-glucose  complex,  and  then  the  FAD 
is  reduced  according  to 

C6Hi206(glucose)  +  FAD  =  CgHioOg  +  FADH^.  (1) 

The  reduced  form  of  the  enzyme  is  fast  oxidized  by  oxygen, 
producing  hydrogen  peroxide,  and  it  restores  the  initial  state 
of  the  enzyme  molecule 

FADH2  H-  O2  =  FAD  H-  H2O2 .  (2) 

In  the  second  substrate  reaction,  the  generated  gluconolac- 

tone  (CgHioOfi)  is  hydrolyzed  spontaneously  to  gluconic  acid 
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FIG.  3.  (Color  online)  (a)  Device  response  as  a  function  of  glucose  concen¬ 
tration,  taken  on  the  second  day.  (b)  Device  response  as  a  function  of  the 
reference  gate  voltage  in  16  mM  glucose  solution,  taken  on  the  second  day. 
(c)  Calibrated  surface  potential  change  vs.  glucose  concentration,  taken  on 
different  days.  All  solutions  contain  50  mM  of  NaCl  and  81.8  mM  of  po¬ 
tassium  ferricyanide. 

CfiHioOe,  +  H2O  =  C6Hii07(glucose)  +  H".  (3) 

Therefore,  one  glucose  molecule  yields  one  hydrogen  ion, 
and  it  changes  the  local  hydrogen  ion  activity  of  the  solution. 
In  our  earlier  work,  we  have  shown  that  the  top-down  fabri¬ 
cated  sensors  are  sensitive  to  pH,  which  quantifies  the  hy¬ 
drogen  ion  concentration.  Changes  in  the  local  hydrogen  ion 
concentration  alter  the  surface  potential  and,  hence,  the  elec¬ 
tric  field  which  modulates  the  conductance  of  the  BioFET. 
This  field  effect  is  further  amplified  by  adding  ferricyanide 
ions  {[Ee(CN)6]3-}: 

EADH2  H-  2[Fe(CN)6]^-  =  FAD  -b  2[Fe(CN)6]^-  -b  2H+, 

(4) 

Fig.  2(b)  confirms  this  effect.  An  81.8  mM  potassium  ferri¬ 
cyanide  solution  has  the  sensitivity  of  16  nS/mM,  while  a 
50  mM  potassium  ferricyanide  solution  gives  the  sensitivity 
of  10  nS/mM.  So  higher  concentration  of  potassium  ferri¬ 
cyanide  yields  higher  sensitivity. 

To  evaluate  the  device  performance  and  degradation  as 
the  function  of  calibrated  surface  potential,  similar  to  the 
macroscopically  large  ISFET  experiments,  we  have  mea¬ 
sured  the  calibrated  plot  of  surface  potential  change.  The 
surface  potential  is  offset  from  the  reference  electrode 
potential  because  of  contact  potential  drop  and  electrolytic 
screening  effects.  However,  the  change  in  surface  potential 
induced  by  oxidation  reaction  is  equal  to  the  change  in 
reference  electrode  voltage  AV^g  applied  to  achieve  the  same 
device  response  in  the  solution.  This  is  consistent  with  our 
previous  experiments,  in  which  we  showed  the  equivalence 
between  changing  reference  electrode  voltage  and  changing 

analyte  concentration  in  solution.*"^  Figure  3  shows  a  plot  of 
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the  calibrated  surface  potential  change  versus  glucose  con¬ 
centration.  First,  we  perform  real-time  differential  conduc¬ 
tance  measurements  by  changing  glucose  concentration  [Fig. 
3(a)].  Then,  at  the  same  V  in  the  16  mM  glucose  solution,  we 
perform  a  rapid  scan  of  to  obtain  a  similar  plot  [Fig. 
3(b)].  The  differential  conductance  of  the  device  is  fitted 
with  an  exponential  growth  function  of  the  reference  gate 
voltage, 

=Co  +  G.»p(^).  (5) 

rg 

and  the  value  of  is  found  to  be  77  mV  for  the  listed 
sample.  This  fitting  parameter  implies  how  the  device  con¬ 
ductance  changes  according  to  surface  potential  change  by 
adding  analyte  concentration.  Lower  value  of  Vt  gives  higher 
sensitivity  of  the  device.  Using  the  data  from  these  two 
measurements,  we  can  plot  the  calibrated  surface  potential 
change  versus  glucose  concentration  by  converting  device 
response  into  an  equivalent  surface  potential  change.  In  this 
step,  we  convert  each  differential  conductance  of  the  device 
in  different  solutions  into  an  equivalent  reference  gate  volt¬ 
age  by  using  the  curve  in  Fig.  3(b),  and  we  count  the  voltage 
difference  between  different  solutions.  The  final  plot  of  cali¬ 
brated  surface  potential  change  as  the  function  of  glucose 
concentration  is  shown  as  the  second  day  curve  in  Fig.  3(c). 
From  the  first  day  curve,  the  slope  of  160  mV/decade  change 
of  glucose  concentration  is  observed.  The  slope  decreases 
after  operation  of  several  days,  implying  degradation  of  de¬ 
vice  performance.  This  degradation  is  attributed  to  a  side 
reaction  of  the  released  hydrogen  peroxide,  which  can  di¬ 
rectly  oxidize  and  deactivate  of  the  enzyme  on  the  device 
surface.  Thus,  the  application  of  the  device  is  limited  to 
short-term  use.  However,  long-term  use  could  be  achieved 

by  inducing  a  secondary  active  electrode  near  nanowires  to 

1  2 

eliminate  degradation. 

Taken  together,  our  studies  show  that  the  nanowire  sen¬ 
sors  provide  an  excellent  platform  for  disposable  glucose 
sensors.  In  order  to  implement  the  sensors  for  long-term  im¬ 
plantable  use,  a  solution  to  the  degradation  problem  is 
needed.  One  suggestion  first  proposed  for  macroscale  ISFET 


using  a  control  electrode  to  counteract  the  effect  of  the  per¬ 
oxide  appears  to  be  promising.  It  would  be  important  to  ex¬ 
tend  the  performance  of  nanowire  sensors  by  exploiting  a 
similar  idea. 

In  conclusion,  we  have  fabricated  silicon  nanowires  as 
glucose  biosensor.  In  the  differential  conductance  measure¬ 
ment,  the  device  response  has  a  linear  range  for  glucose  con¬ 
centration  between  0.5  and  8  mM  upon  initial  use.  While  the 
calibrated  surface  potential  versus  glucose  concentration  plot 
gives  the  maximum  slope  of  160  mV/decade.  The  device 
also  performs  with  a  reasonable  response  time  of  3-5  min. 
The  linear  glucose  detection  range,  fast  response  time,  and 
detection  limit  provide  a  pathway  to  fabricate  high-density, 
high-performance  nanoscale  glucose  biosensors  that  can  be 
integrated  with  silicon-based  circuitry. 
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